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The electric microfield distribution (MFD) at an impurity ion is studied for two-component (TCP) 
electron-ion plasmas using molecular dynamics simulation and theoretical models. The particles are 
treated within classical statistical mechanics using an electron-ion Coulomb potential regularized at 
distances less than the de Broglie length to take into account quantum-difi'raction effects. Corrections 
to the potential-of-mean-force exponential (PMFEX) approximation recently proposed for MFD in 
a strongly coupled TCP [Phys. Rev. E 72, 036403 (2005)] are obtained and discussed. This has 
been done by a generalization of the standard Baranger-Mozer and renormalized cluster expansion 
^ ' techniques originally developed for the one-component plasmas to the TCPs. The results obtained 

^ O I for a neutral point are compared with those from molecular dynamics simulations. It is shown that 

. the corrections do not help to improve the PMFEX approximation for a TCP with low ionic charge 

Z. But starting with Z > 5 the PMFEX model is substantially improved and the agreement with 
numerical simulations is excellent. We have also found that with increasing coupling the PMFEX 
approximation becomes invalid to predict the MFD at a neutral point while its corrected version 
agrees satisfactory with the simulations. 

p ^ PACS numbers: 52.27.Gr, 52.27.Aj, 52.65.Yy, 05.10.-a 
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I. INTRODUCTION 



C/3 

^ ' The determination of the electric microfield distribution at a neutral atom or charged impurity ion (radiator) in 
' ^ ^ a plasma is an important tool for the understanding of many spectroscopic experiments [H, [2j] . Since the pioneering 
J>.^- work of Holtsmark Q , who completely neglected correlations between the particles (ideal plasma) , many efforts have 
■ been concentrated on an improved statistical description of the microfield distribution. The first theory which goes 
Qh. beyond the Holtsmark limit and which is based on a cluster expansion similar to that of Ursell and Mayer [3| was 
' ■ developed by Baranger and Mozer 0, Q . In this approach the microfield distribution is represented as an expansion 
, in terms of correlation functions which has been truncated on the level of the pair correlation. The latter is treated in 
^ ' the Debye-Hiickel form which corresponds to the first order of the expansion in the coupling parameter. The theory by 
' Baranger and Mozer was improved by Hooper 0, @| . They reformulated the expansion of the microfield distribution 
xj in terms of other functions by introducing a free parameter which was adjusted in such a way to arrive at a level 
where the resulting microfield distribution did not depend on the free parameter any more. However, it was argued 
. that such a method is only valid at small coupling parameters, where the correction to the Holtsmark distribution, 
' corresponding to the first term in the series, is small. The first theory capable to provide reliable numerical results 
. for strongly coupled plasmas, known as adjustable-parameter exponential approximation (APEX), was proposed by 
Iglesias and co-workers @, [l^. It involves a noninteracting quasiparticle representation of the electron-screened ions. 
This phenomenological but highly successful approximation is based on a parameterization of the electric microfield 
produced on a radiator with the Debye-Hiickel-type screened interaction with unknown screening length. This free 
parameter is then adjusted in such a way to yield the correct second moment of the microfield distribution. Afterwards 
^ the APEX model was substantially improved for neutral radiators using renormalized correlation functions and electric 
d ' fields ^Tj. (We refer the reader to Ref. for a recent review). 

Most of this work was done on one-component plasmas (OCP) and thus neglects the influence of the attractive 
interactions between electrons and ions. In this paper we study the microfield distribution in a two-component 
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plasmas (TCP). This has been done previously in Refs. [l^, [l^ for partially degenerate electrons. In Ref. [l^ 
the low-frequency component of the microfield was calculated within the linear response treatment taking strong 
correlations into account via local field corrections. Also the problem of attractive interaction has been considered 
for a single but highly charged impurity ion immersed in an electronic OCP [15]. In the recent papers [l^, [l3| 
we already studied strongly coupled systems, i.e., a highly charged radiator in a TCP of classical (nondegenerated) 
and strongly correlated particles beyond a perturbative treatment. As in Ref. the theoretical scheme (PMFEX) 
presented in Refs. [IB, Il3| is based on the potential-of-mean-force approximation. It exactly satisfies the sum-rule 
requirement arising from the second moment of the microfield distribution without introducing adjustable parameters. 
Comparisons of the PMFEX calculations with the molecular-dynamics simulations for the electric fields at the highly 
charged points relevant for laser-produced plasmas show, in general, an excellent agreement even for large coupling 
p^ . However, the results of the PMFEX for neutral points have not yet been comprehensively tested. Here, the 
PMFEX scheme may be less accurate , similar as it occurs in the APEX where further improvement for microfield 



distribution at neutral point is needed, see Ref. 11|. Moreover, in a regime dominated by small local fields and hence 
by small local electronic density the PMFEX deviates from the molecular-dynamics simulations also for charged 
radiator This feature has been clearly observed for a single ion embedded in an electronic OCP in Ref. [l5l |. 

Therefore, it is of interest to improve the PMFEX model and provide a tool for calculating the corrections to it 
when required. Our purpose here is to introduce the key ideas of PMFEX in the standard Baranger-Mozer cluster 
representation for the microfield distributions [B, @| . In this way some contact between PMFEX and standard small 
plasma-parameter theories is established, while also providing PMFEX as the leading term in a series from which 
corrections can be calculated explicitly. Another important ingredient of our treatment is the electron-ion attractive 
interaction which drastically changes the physical properties of the system as compared to classical OCPs (see, e.g., 
Ref. (15|]). But the thermodynamic stability of a TCP requires at least some quantum features for the electron- ion 
interaction at short distances. This required minimum of quantum features is here taken into account by using a 
regularized ion-electron potential [H, [3 113 7 which then enables the application of classical statistical mechanics and 
classical MD simulations. 

The paper is organized as follows. In Sec. HIl we define the basic parameters of interest for a TCP. The theoretical 
schemes applied previously to either electronic or ionic OCPs are generalized to TCPs in Sees. IIIII and IIVI In 
particular, starting with the canonical ensemble formulation of the microfield distribution in the TCP the Baranger- 
Mozer treatment is reviewed and the structure of the resulting cluster series is discussed in Sec. IIIII Next, the basic 
assumptions of the PMFEX are briefly given and motivated in Sec. IIVI Here a formal relationship of the PMFEX to 
the Baranger-Mozer series is established and the first two terms of a renormalized cluster series are given explicitly. 
Furthermore, to test obtained theoretical results in Sec. |V] the correction terms are calculated at a neutral point and 
compared with the results from classical molecular dynamics simulations. The results are summarized in Sec. I VII 
The exact second moment of the microfield distribution at a neutral point is considered in the Appendix El 

II. BASIC PARAMETERS FOR THE TCP 

We consider a neutral and isotropic two-component electron-ion plasma consisting of ions and electrons at a tem- 
perature T in a volume fi. The particles are assumed to be classical and pointlike. The average densities, charges, 
and masses of the ions and electrons are n^, rig, and Ze, — e and mi, m, respectively. We assume that the density of 
radiator ions or atoms are small, hr <§; rii-e and thus consider only one radiator ion with charge Zne in our calculations 
(throughout this paper the index R refers to the radiators). Because of the charge neutrality we have rig = riiZ. 

We now introduce the Coulomb coupling parameters Tap- Introducing the Wigner-Seitz radii, i.e., the mean 
electron-electron, electron-ion, and ion-ion distances through the relations, a^^ — Anrie/S, a^^ — Ann/S, and a^^ — 
Aimi/S (where n = + rii is the plasma total density) these parameters are defined as 

r r T (1) 
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respectively, where e|. = e^/Aneo. Note that 

J- ee — 1 /rt 1 ^ a — 1 /q • \^ ) 

[Z-^iZ+l]]^^^ [Z + Xf'^ 

In a hydrogen plasma with Z — \ we obtain Fee = Va = 2~^I'^V while in a plasma with highly charged ions (Z 3> 1) 
Vii = Z^/^Tei ^ Fee — ^ei/Z. So, in the TCP with uncorrelated electrons the ions may be strongly correlated. 

The Holtsmark field Eh for a TCP is given by E^J'^ = E^J^ + E^J^ (see Ref. ^ for details), where Ehc and Eni 

are the electronic and ionic Holtsmark fields, respectively, Ehc — Cep /a\, Ehi — CZep/af with C = {8tt/25)^^^ and 



ep — e/ineo- Since Ehc = Z^^/'^Ehi the electronic and ionic components of a hydrogen TCP contribute equaUy to 
the Holtsmark field. For a completely ionized TCP with highly charged ions the ions dominate Eh- The definition of 
the Holtsmark field Eh for a TCP is equivalent to the obvious relation n = rif, -\- rii and can be represented as 
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with an effective charge Z. 

For the regularized pair interaction potential e^gqaqpUap (r) with a, /3 = e, i, i?, qe — —1, ft = qn — Zfi, we here 
employ 



(4) 
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The cutoff parameters 5aj3 are related to the thermal de Broglie wavelengths, Sap = / IJ-apknT^ , where ij^ap 
is the reduced mass of the particles a and (3. For large distances r > Sa/i the potential becomes Coulomb, while 
for r < Sap the Coulomb singularity is removed and Uap{0) = l/daf3- By this the short range effects based on the 
uncertainty principle are included [H, [l^, [131 • We point out that the dependence on the physical parameters density 
and temperature is implicitly contained in the parameters Tap and 6a/3- 



III. GENERALIZATION OF THE BARANGER-MOZER FORMULATION TO THE TCP 

A. Preliminaries 



The electric microfield distribution (MFD) Q (E) is defined as the probability density of finding a field E at a 
charge Z^e, located at Tq, in a TCP with Ni ions and Ng electrons. This system is described by classical statistical 
mechanics in a canonical ensemble of (Nj + Ng + 1) particles, and temperature T. The normalized probability density 
of the microfield E is then given by [l^, [13] 



Q{E)^^ f e-^-^(^-^-"-°)(5(E-E(Te,7-,ro))drodred7-, 



(5) 



where (3t — ^/ksT, and 7^ = {ri, r2...rAr^}, % — {Ri, R2...RAr. } are the coordinates of electrons and ions, respec- 
tively, and 
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is the canonical partition function. U (7^,7i,ro) is the potential energy of the configuration given here by 
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in terms of the pair interaction potentials Uq/3 (r) and UaB.{r), where a, f3 = e,i, ri*^-* = Tq, Va' = Ra. In Eq. ([7]) 
the first sum is restricted to a 7^ 6 for identical particles, a = p. The total electrical field E (7^, 7^, rg) acting on the 
radiator is given by the superposition of electronic and ionic single-particle fields 



E (Te, 7-, ro) = -J^^oU = Yl ^^(^0 ~ ^'') + E''(^o - Ra) 



(8) 



As Eq (r) = ^Ea{r), we obtain for the electronic and ionic single-particle fields Ef.{r) = eF^efl (»'), Ei{r) = 
—Zepu'^jj (r), where the prime indicates derivative with respect to r. 

The spherical symmetric interaction between plasma particles and the isotropy of the system allows to introduce 
the normalized microfield distribution P{E) = AttE^Q{E). It can be reexpressed in terms of the Fourier transform of 
Q{E) through 
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(9) 



and 



(10) 



where (...) denotes a statistical average and ^0(2;) = sinx/x is the spherical Bessel function of order zero. The 
coefficients of the expansion of the function T{k) at k — > yield the even moments of the microfield distribution, 
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The similar expansion for the function C (k) defined by T(k) 



yields 



(11) 
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Therefore the Fourier transform of the MFD can be interpreted as a generating function for microfield even moments. 
Equations ((5|)- (fT^ describe the total MFD at the position rg of the radiator generated by both the statistically 
distributed ions and electrons of the TCP. Since we are interested to calculate the MFD, Eq. in an infinite system 
the statistical average of any quantity becomes translationally invariant with respect to Tq, and the location of the 
test charge may be taken as the origin without loss of generality. 



B. Baranger-Mozer cluster expansion for the TCP 

In this section we generalize the Baranger-Mozer (BM) cluster expansion technique originally developed for one- 
component plasmas (see, e.g., [1, H, 0, @l) to the classical TCPs. This method results from two transformations of 
Eq. ^ . The first is motivated by the fact that the required average in Eq. (fTU|) is the product of electronic and ionic 
single-particle functions, e*'^ -^^^""") and e"''^*'^"', which have a value close to one over some volume of the system. 
This suggest a first transformation to the set of single-particle functions 

xi")(/.) = e'''-=^°("-»°')-l. (13) 

Using this transformation the exponential factor in Eq. (|10p becomes (for simplicity we drop the coordinate Tq of the 
radiator) 



exp [^K ■ E {%,%)] - n [1 + ('^)] n [1 + ('^)] • (14) 
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The introduced single-particle functions Xa have the advantage to be zero over some volume. The spirit of this 
transformation to the functions Xa^^ is similar to the use of Mayer's /-functions for thermodynamic properties of 
gases |4j]. Substituting Eq. (jl4p into Eq. PH)) leads in the thermodynamic limit directly to the series 

r(K) = 1 + ^5:^ / x^^\K)xtH>^)-x^aH^)Qi''Kr^-))dT^-^ (15) 
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where rj^^ = {r|"\ ...ri^^} and dTj") = dv^{'Ur'^^\..dv^^\ Here ^^"^(ri"^) and ^l^^ ^/'^ ) are equilibrium 
correlation functions. The first quantity represents the probability density for a particles from plasma species a at 
rj"-*, r2"\...ri"\ and the test particle at the origin. The second one describes the correlations between a electrons at 



ri, V2,---Va, and b ions at Ri, R2,...Rb involving the test particle. The range of integration for each term in Eq. (fli 
now restricted by the xi"' (K)-functions. However, this restriction is not uniform with respect to and particularly 

for large values of k the functions xi"' ('^) can differ from zero over a correspondingly large volume. Consequently, a 
second transformation is desirable, T{k) — e~^'-'^\ where by a standard theorem of equilibrium statistical mechanics 
(see, e.g., [2l|), C{n) is determined from Eq. (fT5|) as 

00 ^ °° a °° b 

^ (-) = - E E ^'^i"^ («) - E ^ E ^f^' (-) • (16) 

a a—1 a— 1 b—1 



Here ha and /i^^*' (k) are given by 
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and ^i"'' and l^^^^ are the Ursell cluster functions for TCP associated with the set of usual correlation functions 
^a"^ and Q^ab^ , respectively. The functions fi"^ for electron-electron and ion-ion interactions are expressed by the 
correlation functions ^/a"'' in the same manner as in the case of corresponding OCP (see, e.g., @, H, 0|)- Therefore, as 
an example we consider here explicitly only the functions i^^^^ which involve electron-ion interactions, 

4f (ri,Ri) - (ri,Ri) - 5^^' (ri)^« (Ri) , (19) 

4'^(ri,Ri,R2) = (ri,Ri,R2)-ei'^^(ri)gf(Ri,R2)-^f^(Ri)gif (ri,R2) (20) 

-ef ^ (R2) e^r^ (ri, Ri) + 20^^^ (ri) e« (Ri) ^ (R2) , 

4f (ri, r2, Ri) = CJ^r^ (ri, r2, Ri) - (Ri) G^^^ (ri, r^) - g['^ (n) g^f^ (^2, Ri) (21) 

(r2) air^ (ri,Ri) + 2^1^^ (ri) 0^^) (rs) (Ri) 

etc., for pair and three-body correlations, respectively. The higher order Ursell functions involving ^ 4 particles 
can be constructed similarly. All single-particle Ursell functions are here expressed by the ordinary pair correlation 
function between radiator and plasma particles of species a, i'f^ (r) = ^("^ (r) = gaR{r) [lB|- The significant 
difference between Q and i correlation functions is that the cluster functions £ vanish when any members of the 
particles are sufficiently far apart, whereas the ordinary correlation functions Q have not this property. In particular, 
for completely uncorrelated system of particles the many-body functions Q tend to unity while ^ — > beginning with 
i'fi^ and Hence, the range of the integrals in Eq. ([TC]) is controlled by both the functions xi"' (k) and the Ursell 
functions. The later restricts the integration to volumes characterized by the correlation length which depends on 
the thermodynamic-state condition but is independent of k. Qualitatively, therefore, the BM formalism provides a 
series representation whose terms are controlled by the range of xi""* (k.) for small k and by the range of the Ursell 
functions for large k. 

Equations (fT6|l -(fT8 |) together with the relation T(/«) = e"'^''^-' constitutes the generalization of the Baranger-Mozer 
cluster expansion technique to the TCPs. In contrast to the BM theory developed for OCP Eq. now involves 
terms which are responsible for electron-electron, ion-ion and electron-ion attractive interactions. The BM result is 
easily recovered from Eq. (fTH]) by neglecting the last term as well as one of two terms responsible for the interactions 
between identical particles. In addition, for ideal TCPs (Holtsmark limit) the second sum in Eq. (fT6|) and all terms 
with a ^ 2 in the first sum vanish, and only the term with a = 1 in the first sum contributes to the MFD. It is easy 
to see that this term coincides with the result of Ref. [16] derived for ideal TCP. 

From a practical point of view, it is too difficult to calculate correlation functions of higher order than ^^2""* and ^j;^*"* . 
For weakly coupled plasmas, the Ursell functions of order A -I- 1 are typically of the order of the plasma parameters to 
the power A. Therefore Eq. is usually truncated at first, second or in some cases at third order. If the correlation 
functions are expanded as well with respect to the plasma parameters and truncated at the corresponding order the 
MFD with Eqs. ((IHll-dlHl) can be evaluated analytically (see, e.g., Refs. |E H, 0, H, HI, HI] ) . Obviously, such a method 
may be successful only for weakly coupled plasmas. A possible extension to the strongly coupled regimes requires 
that the BM representation sufhciently rapidly converges to allow for a truncation after the first two terms. This 
suggests, for instance, that the BM series could be improved if the bare single-particle electric fields in the x-functions 
are replaced by screened fields representing the effects of strong correlations in a plasma. The formal procedure for 
carrying out such a " renormalization" has been developed in Ref. [2^ and has been previously employed for OCPs 
[l]| . In the next Sec. IIVI we extend this renormalization procedure to the two-component plasmas. 



IV. RENORMALIZATION OF THE BARANGER-MOZER CLUSTER SERIES 



A. PMFEX as independent-particle model 



In Refs. [ia, we introduced the potential of mean force exponential approximation (PMFEX) which links the 
MFD to the pair correlation functions. To make further progress we briefly outline here the basic concepts of the 
PMFEX approximation. We show that similarly to the APEX (see, e.g., H, [13, |Tl[) PMFEX is also an effective 
independent-particle model. Based on this assumption the MFD within PMFEX is derived quite simply. If all 
interactions between the plasma particles, except those with the radiator, are neglected, then the Ursell functions fi"^ 
vanish for a > 2 and i'"^^'' {Ta'^^ , 7^^*^ ) = for arbitrary a and b. The BM series (fTB)) then reduces to only the leading 
term, 

in) = -J2 I x["^ a^R (ri) dri. (22) 

The superscript (0) denotes the corresponding quantity without interactions among plasma particles. The PMFEX 
approach retains the independent-particle form of (I22p . 

{k) ^C{K)^-J2na f 4"^ («) 9*aR (ri) dr„ (23) 

with the assumption that the important effects of correlations can be accounted for by an effective pair distribution 
function g*^ (ri) and a screened field £a (ri) replacing the single-particle field Eq (ri) in xf^ ('*)j 

V-i"^ (k,) - £"'■^-('■1) - 1. (24) 

Some constraint must be imposed to determine 5*^ (ri). As in the case of APEX [il this is a requirement that the 
effective " quasiparticle" field due to the effective charge density at ri is equal to the corresponding exact field 

9aR (i-i) £a (ri) = g»R (ri) E„ (n) . (25) 

Assuming the spherical symmetric interactions with £a (r) — -£a (f) and (r) = -Eq, (r), this condition then define 
the PMFEX model as 
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Jo 



- jo {nSa (r)) 



g^R (r) r'dr. (26) 



In practice gaR (r) is calculated from the hypernetted-chain (HNC) integral equation [2^, [2^ extended to the TCP 
(see, e.g., Refs. [laflTj). 

In the APEX model originally developed for the classical ionic OCP with bare Coulomb interaction to determine 
the effective field a second requirement is needed. In Rcf. [9] the effective field £{r) is assumed as a Debye-Hiickcl-type 
screened interaction with unknown screening length. This free parameter is then adjusted in such a way to satisfy 
the exact second moment requirement. In contrast to the APEX the effective fields £a (r) in the PMFEX model are 
obtained automatically (see [16| for details) employing the thermodynamic perturbation theory [23|- We recall here 
these effective fields for completeness. For simplicity assuming Coulomb potential for the repulsive interactions, i.e. 
^ee — Sii — (clectron-electron and ion- ion interactions) and regularized one (see Eq. Q) for electron-ion interactions 
these fields read [11] 

£. (r) = (r) ( I + r (p) - g,, (p)] p'dp] , (27) 
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where Ea (r) = qaep/r"^ are the bare single-particle fields. Based on the derivations given in Ref. it is straight- 
forward to show that the effective fields £a (f ) can be alternatively represented as the logarithmic derivative of the 
radial distribution functions (RDF) gaR (r) as 

£.ir) = ^^[\ngaRir)] (28) 
Zre or 



which is known as the potential of mean force (PMF) approximation proposed by Yan and Ichimaru (see also 
[25II). Using Eq. (fT^ we expand around a value k — and find the second moment of the MFD in the PMFEX 
model 

/>oo 

(^')PMFEX = E 4^^"" / (^) ^''^^^ (29) 

Then, introducing Eq. (^5]) in Eq. (^5]) automatically satisfies the exact sum rule, (£'^}pmfex — {E^), without any 
adjustable parameter. Here (i?^) is the exact second moment of the MFD derived in Ref. [l^ for the TCP, see 
Eq. (IA1|) . Thus, if the gaR{r) arc known the MFD within the PMFEX model and with the exact second moment can 
be calculated using Eqs. ([26]) . and ([M]) . 

The PMFEX effective field, Eq. is valid for charged radiator with Zji ^ 0. For a neutral radiator, i.e., in 

the limit Zn — > 0, the PMF ansatz (|28p is not applicable, as in this case the RDFs tend to unity, gaR 1, and 
^T^QaRif) 0. Based on the derivations given in Refs. 16, 17j (see also Appendix [X] for details) the correct limit 
Zr can, however, be done and results in the effective field of Eq. (|27|) . by setting gaR (r) = 1. 



B. Renormalized cluster series 



The comments at the end of Sec. IIIIBI and the success of PMFEX indicate that the Baranger-Mozer series could 
be improved if the single-particle fields in the functions xi"'' (*«) is replaced by a screened fields representing the 
effects of correlations in strongly coupled plasmas. In this spirit, a new functional series is obtained in terms of the 
renormalized functions ■(/'i"'' (k) of Eq. by the definition £* [t/jq; k] = C [xa', 1^] which is easily obtained from the 
BM series Eqs. (|16p - p^ and the functional relationship between xi"-* and Introducing the functions TZa (r) 

by the relation £a (t) — Ti-a {t) Ea (r) (we assume spherical symmetric interactions) with the effective £a (t) and 
single-particle Ea (r) fields, the functional relation between xi""* and V'i"^ is established as 
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To obtain the new function C* [ipa'iK] which must be multi-linear with respect to ipa (cf. Eqs. (|16 p -p8 p ) we expand 
the Xa"^-functions in Eq. ([501 with respect to i/ii^^-functions 
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Then elimination of xi"'' on the right side of equation C* [ijja] k] = C [xa] 1^] using Eq. (|3ip yields the desired renor- 
malized cluster series, 



[^.; = - E E ^^^^ ('^) - E ^ E ^Hif (.) . 
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Here 



(33) 



H. 



(ei) 
ab 



{k) = / (k) (k) ...v^i^' (k) (k) (k) ...4:^ (k) 



(34) 



'(0 



with the generahzed UrseU functions L^"'' and which are recognized as the functional derivatives of H^"^ (k) and 
^ab ^ (k), respectively. Explicitly the first and second order generalized Ursell functions are found to be 

^(.) (35) 
TZa (ri) 
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T^a (ri) 7^a (r2) ria (ri) 
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which result the first three terms in Eq. 
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The first term of Eq. is seen to be precisely PMFEX (cf. Eq. and ^^). The factor Tla (r) occurs automat- 
ically here from the renormalization and eliminates the assumption ([25|) of PMFEX. 

Finally the angular integration in Eq. (|38p can be performed using spherical harmonic expansion [28| . This yields 

£ (k) = ^PMFEX (k) + AZ: (k) , (39) 



where the first term is the PMFEX result, Eq. ([26| . and the second term represents the corrections due to the 
renormalization, 
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Here Cq (k) is the /Ipmfex (k) given by Eq. ([26]) but with ga/j (r) = 1, 

In Eq. (j40p the term 5a/3 (^) (with = 6,2) is the static structure factor for the two-component plasma, 

Sa0 (k) = Sai3 + 'inuaf} [gafj (r) - 1] jo {kr) r'^ dr (42) 



(41) 



with Saa — 1, Sei — 0, — 7 ^ei = u = 71^. -\- and Gap (/^, k) is defined by 

00 

G^f, {k, k)^Y. (-1)' + 1) ■t'' k) {k, k) , (43) 



1=0 



ji {kr) [ji {K£a{r)) - 5io\ , 

'^a V J 



(44) 



TABLE I: The critical values ac{Z,5) for some values of the ion charge Z and three values of 5 — 0.1a, 0.2a and 0.4a. 



s 


TT + 


T -3 + 




N' + 


A 1 13 + 


Ca ^ 


0.1a 


6.87 


7.59 


7.67 


7.62 


7.59 


7.63 


0.2a 


6.55 


7.68 


7.52 


7.13 


6.66 


6.48 


0.4a 


8.90 


12.73 


10.76 


9.39 


6.75 


6.00 



Also, ji (x) is the spherical Bessel function of order I. It should be noted that for deriving the renormalized series 
Eq. ((5^ we do not use explicit functional form of £a{r). Hence, in spite of its important role in the theory, the precise 
functional form of the effective field Sair) may remain arbitrary. Moreover, it is straightforward to show that the 
generating function Eq. as well as its renormalized version Eq. ([5^ yield the exact even moments of the MFD 
according to Eq. (|12p independently of the functional form of £a{r). In the present context it might appear more 
reasonable to choose ^^(r) to improve convergence of the renormalized series Eq. p2p . A similar procedure is used 
in thermodynamic perturbation theory where the corresponding parameters of the leading term are chosen to make 
the next order terms vanish . But this is not possible without making £a (r) a function of k (see similar discussion 
in Ref. [ll'| in the context of APEX). Here we consider for £a{r) the mean- force field as derived in Ref. for the 
TCP (see also Eq. ([28l) as well as Eq. ([27l) for the repulsive Coulomb interactions) and which appears automatically 
from the second moment condition d^£* /dK^\^_^ — \{E'^) for the renormalized function C* (15^ . The theoretical 
scheme resulting from Eqs. is abbreviated as PMFEX-I-. It agrees for neutral points quite well with the 

molecular dynamic (MD) simulation results, as we will show in the next section. It is straightforward to show that 
the correction A£(k) in Eq. ([5^ behaves as ~ k'' at k -+ and, therefore, does not contribute to the second moment. 
Then the quantity receives contribution only from the PMFEX term in Eq. ([39]) . 

V. APPLICATION TO THE NEUTRAL RADIATOR 

In the case of neutral radiator, Z/j — 0, the plasma-radiator correlation function is gaR (r) = 1 and Eq. (|40p is 
simplified. In particular, the last term vanishes since >Cpmfex (k) — Cq (k) for (Jqa (r) = 1. Thus, the generating 
function C (k) of the MFD in the PMFEX+ model is given by C (k) = £o («) + A£ (k), where Cq (k) is the ordinary 
PMFEX expression, with gaR{r) = 1, see Eq. (|41]) . The higher order corrections are involved in the second term, 
AC (k). The functions TZa (r) are obtained from Eq. (P7)) . 

TZa (r) = 1 + inUa / [gaa{p) " gei{p)] P^dp (45) 

Jo 

which in the limit of small plasma-parameters [T^e — > 0) can be evaluated explicitly using the Debye-Hiickel approxi- 
mation. Tie {r) = Tli {r) ~ ( j -I- 1) e^^'l^. Here = e^ksT / {Zne^) is the Debye screening length. Thus similar to the 
standard PMFEX, the PMFEX+ approximation links the MFD to the RDFs. To obtain exphcit results for the MFD 
the corresponding RDFs and the static structure factors must be determined first. This is done by solving numerically 
the hypernetted chain (HNC) integral equations [23, HH] for the TCPs under consideration (see Ref. [ifif and references 
therein for more details). For simplicity we here assume bare Coulomb interaction for electron-electron and ion- ion, 
i.e. See — and da ~ 0, and a regularized ion-electron interaction with a parameter 5ei = 5 fixed to 5 ^ 0.2a or 
5 = 0.4a, where a = (47rrt/3) is the Wigner-Seitz radius. For this kind of TCPs, the HNC method has already been 
extensively tested and evaluated by comparison of the resulting RDFs with those obtained by classical MD simulations 
p^ . For the numerical solution of the HNC scheme the dimensionless parameter Oei — Ze^Uei (0) /ksT = Tei (a/ 5), 
i.e. the maximum value of the electron-ion interaction energy in units of fesT, plays an important role (see also 
0, [iB|)- Within our numerical treatment of the HNC equations a parameter regime with Oei < cFciZ, 5) is accessible, 
where the critical values CTc for 5 = 0.1a, 0.2a and 0.4a and the different studied TCPs (H+, Li^+, 6^+, N'^+, k\^^+ 
and Ca^*^^) are given in Table HI Beyond this value the HNC numerical procedure does either not converge or ends up 
in unphysical solutions. A similar behavior has been reported in Ref. [15] for the case of an ion embedded in electrons. 
With the RDFs and Sap (k) provided by the HNC scheme the MFD P{E) in the PMFEX4^model is then calculated 
via Eqs. ^ and ([39|) - ([45|) by standard numerical differentiation and integration methods 30]. In practice it appears 
sufficient to terminate the sum in Eq. (|43p at ? = 5. 

To test the PMFEX-f approximation and improvements compared to the standard PMFEX approach the calculated 
MFDs are confronted with those obtained by classical MD simulations. In the MD simulations the classical equations 




FIG. 1: RDFs 5a/3(r) for (left panel) and Ca^°+ (right panel) plasmas with 5 = 0.2o, Tee = 0.25 and Tee = 0.03, 

respectively. The lines correspond to the HNC approximation while the symbols denote the MD simulations. The different 
lines represent gee (green), gu (red) and gei (blue). 




r/a 



FIG. 2: Same as Fig.[T]for N'^+ plasma with <5 = 0.2o, Fee = 0.15 (left panel) and Tee = 0.18 (right panel). 



of motion arc numerically integrated for Ni ions and Ne — ZNi electrons interacting via the regularized potentials 
Uafsif)- Such MD simulations have already been extensively tested and successfully applied for investigations of the 
dynamic properties [Sll, [H, [H, [H, [1^ Ha] and the MFDs [H, |T3] of a TCP with regularized potentials (see also 
Rcf. [38] and references therein). For our present considerations we investigated the specific cases of H"'", Li^"*", B^+, 
, Al^'^^ and Ca^°+ TCPs with symmetric (hydrogen) and strongly asymmetric densities of the plasma species and 
certain values of the coupling parameters Taf3 with a, /3 = e, i. These coupling parameters Tap, at a given parameter 5, 
are chosen to avoid the mentioned numerical difficulties with the HNC scheme and the formation of unphysical bound 
states in the classical MD simulations. This implies in particular that for highly charged plasma ions with Z ^ 1 
plasma states with strongly correlated ions {Tu ^1) and strong electron-ion interactions {T^i ~ 1) are accessible, 
while the parameter Fee has to remain quite small, T^e — ^ei/Z <C 1. But in the case of highly charged ions even 
for Fee <C 1 the electrons may be strongly correlated due to nonlinear effects as will be discussed below. The MD 
simulations providing the present MFDs and RDFs have been performed using Nmd = Ni + = 5376 particles. 
Further details on the applied MD technique can be found in Refs. [sj, [H, HI, Hj] . As the theoretical models and the 
numerical solutions depend directly on the coupling parameters Tafs and the regularization parameter i5 we discuss, 
as e.g. in Refs. flEl. fl6l. fl7l|. our results in terms of these parameters rather than in the underlying physical values of 
density and temperature. 

Before discussing the MFD itself, we first consider briefly the RDFs and the validity of the HNC approximation 
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FIG. 3: Normalized electric microfield distributions for a hydrogen plasma with S = 0.4a, Fee = F;; = 1, Cei = 3.15 (top) and 
Fee = Fii = 2.5, (Jei = 7.87 (bottom) as a function of the electric field in units of Eh- The dashed and the solid curves are 
the results of the PMFEX and PMFEX+, respectively. The filled circles represent the MFD from the MD simulations. The 
Holtsmark and standard Baranger-Mozer distributions are also shown as open circles and triangles, respectively. 

for the parameters at hand. For small and moderate coupling where the parameter a^i is well below the critical 
value (Tc, the RDFs as either calculated by the HNC scheme or extracted from the MD simulations are always in 
perfect agreement (see e.g. the corresponding examples given in flG^). Some typical examples for the RFDs at rather 
large coupling, where some deviations between HNC and MD may show up, are given in Figs. \T\ and [H Due to the 
regularization of the ion-electron interaction the RDF geii^) is finite in the limit r — > (not visible in Figs. [1] and 
H]) and can be approximated as gei{0) ~ exp{(Tei/'JZ) with TZ — I + {S/a){3Teiy^^ , see Refs. [l^ [l^. Thus the RDF 
geiir) shows the expected growth of correlations with increased coupling and decreased regularization parameter, 
where the HNC scheme tends to overestimate the value of 5ei(0) (not visible in the figures). For example, for N^+ 
ions g™^(0) ~ 57.4, 5^°(0) ~ 48.8 and 5™^(0) ~ 107.3, g^°(0) ~ 96.4 with ,5 = 0.2a, Fee = 0.15 and Fee = 0.18, 
respectively. At strong coupling deviations also occur in the electron-electron RDF gee (r) at small r (as in the case of 
and N^+). The strong electron-ion interaction increases the electron density around the highly charged ion which 
introduces additional correlations between the electrons. This increases the probability of close electronic distances 
and results in the maxima in geeif) at distances r < a. This effect is obviously again overestimated in the HNC 
approach. 

In Figs. [MHl we next compare the MFDs calculated from the PMFEX (dashed Hues) and PMFEX+ (sohd lines) 
schemes as well as from MD simulations (filled circles) where the electric microfields are scaled in units of the Holtsmark 
field Eh [see Eq. ([3])]. The open circles are the Holtsmark MFDs for a TCP with Coulomb potential and the MFDs 
predicted by the first two terms of the standard BM series Eqs. (fT5| - (PT|) are also shown as triangles. We here focus 
on cases of strong coupling where the PMFEX and PMFEX-t- results significantly differ, i.e. where the higher order 
corrections AC (k) ((4T]) substantially contribute to the generating function £ {k) = Cq (k) + AC (k). As can be seen 
in all the presented cases the MFDs obtained both from the PMFEX and the PMFEX+ treatment strongly differ 
from the standard BM electric field distributions. 

Concerning the agreement with the MD simulations the PMFEX-I- scheme turns out to improve the PMFEX model 
substantially for large charge states of the plasma ions. For H+, Fig. [3l no improvement can be found, and the 
agreement with the MD simulations remains rather unsatisfactory, although the PMFEX(+) methods yield a much 




better description of the MFD than the standard BM approach. For Li^+ and B^+, Figs.[4]and[5l the PMFEX+ scheme 
much better agrees with the MD than the PMFEX, but obviously still deserves further improvement. However, for ion 
charges Z > 5, Figs.[B][51 the MFDs obtained from MD simulations are in excellent agreement with those predicted by 
PMFEX+ (except for some smaller values near the maxima) . The rather large deviation of PMFEX from PMFEX+ 
towards a higher probability of low fiels (in particular in Figs. IH] and [5]) indicates that the PMFEX applied to the 
field at a neutral point does not correctly account for perturber-perturber correlations, and, therefore, underestimates 
the MFDs at large fields. As expected for = all distributions shown in Figs. OH] merge at large electric fields, 
ij = E/Eh ^ 6 with the Holtsmark distribution with the asymptotic behavior Ph {r}) — 1.4QQrj~^-^ . This indicates 
that due to electron-electron and ion-ion Coulomb interactions with 5ee — Su ~ the second moments of the MFDs 
shown in Figs. [3][5]do not exist, see e.g. Eq. (|A7p derived in Appendix \K\ 

To gain some more insight of the features and the range of validity of the PMFEX and PMFEX-f models we 
consider the extreme regimes of Figs. [9] and [10] with further increased coupling parameters. In the PMFEX scheme 
and for large values of Fee the effective fields at large distances behave as Eair) ~ ^e""'" cos(6r). Here a and h are 
some parameters increasing with Fee, whereby the involved screening length does not necessarily coincide with 
the Debye length A. The oscillatory nature of the effective fields drastically changes the properties of the generating 
function Co{k) in the PMFEX approximation (see Eq. (|4ip ). which becomes negative at kEh ^ 1. This is not the 




FIG. 6: Same as Fig.[3]for N'^+ plasma with Fee = 0.1, Fa = 2.56, cr^i = 3.66 and 5 = 0.2a. 




FIG. 7: Same as Fig. H for AI^='+ plasma with Fee = 0.01, Fa = 0.72, aei = 0.67 and S = 0.2a. 

case in APEX (as well as in PMFEX for charged impurity ion, i.e. for the MFD at charged point) where the effective 
field is treated in the Debye-Hiickel form and thus is positive and decreases monotonically with r. Therefore the 
K-integration in Eq. (O diverges and hence the PMFEX model becomes invalid for neutral radiators {Zn = 0) (and 
no curves for the PMFEX can be given in Figs. 151 and fT(7|. But even in this extreme regime the PMFEX+ model 
remains valid and the agreement with MD simulations is still quite good for the lower charge state Z — 7 shown in 
Fig. [H] (although not as good as in Fig. [5]), and almost perfect for the highly charged ions as shown in Fig. [TOl 

As discussed above some deviations between HNC scheme and MD simulations occur in the electron-electron RDF 
geei^) at small distances r < a for and N''+ TCPs, see the left panel of Fig. [T]and Fig. [2l Similarly the static 
structure factor See{k) obtained from HNC deviates from MD data at fca > 1. This may be critical for the accuracy 
of the electronic effective field, Eq. and hence for the generating function C{k) Eqs. (|39 p - ((lT|) calculated from the 
HNC approximation. Since the generating function is involved in the exponential factor e^^^'^^ even small deviations 
in C{k) results some appreciable deviations in the MFD. Apparently the PMFEX+ model could be improved in the 
case of 6^+ and N'^+ TCPs in Figs. E] and [H respectively, using the MD data for the RDF gee{r). 




FIG. 8: Same as Fig. |3]for Csi^''+ plasma with Fee = 0.02, F.t = 2.95, ct„ = 2.03 and 5 = 0.2a. 
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FIG. 9: Normalized electric microfield distributions for N''^ plasma with Fee = 0.15 {Ta — 3.85, a^i = 5.49), Fee = 0.18 
(Fii = 4.61, (7ei = 6.59), and 5 = 0.2a as a function of the electric field in units of Eh- The solid and dashed curves are the 
results of the PMFEX+ model. The filled circles and squares represent the MFDs from the MD simulations. 

VI. DISCUSSION AND CONCLUSION 

In this paper we investigate the microfield distributions in a two-component plasmas with attractive electron-ion 
interactions. Based on the renormalization of the standard Baranger-Mozer cluster expansion technique our objective 
here was to derive the corrections to the PMFEX model proposed recently in Ref. The outlined theoretical 

method is abbreviated as PMFEX+. In general the MFDs predicted by PMFEX for the case of charged impurity 
ion are very accurate and its corrections may be quite small, except the case of strongly coupled hydrogen at small 
regularization parameter 6 As discussed in Ref. the PMFEX is less accurate in the case of the neutral 

point and an improvement of the model is required. We thus focused on testing the predictions of the PMFEX-I- 
approximation for the neutral point based on the HNC treatment of static correlations by confronting it with the 
MFDs obtained from MD simulations. One of the basic assumptions made for all these models is the regularization 
of the attractive electron-ion interaction at short distances to introduce quantum diffraction effects in the employed 
classical approach. For the repulsive electron-electron and ion-ion interactions we assumed a bare Coulomb potential, 
for simplicity. 




As examples the cases of H+, Li'^+, B^+, N^+, Al^^+ and Ca^°"'" two-component plasmas (TCPs) with symmetric 
and largely asymmetric densities of the plasma species were considered. Our treatment is limited to a parameter 
regime with a < <Jc, where the critical values CTc for S = 0.1a, 0.2a and 0.4a are shown in Table HI Within this 
parameter regime the gapij") from the HNC equations agree well with the MD simulations (except the correlation 
function gee{r) for the cases given in Figs. [T] and [5]). Beyond these critical a the HNC equations do either not converge 
or end up in unphysical solutions. A further increase of the coupling parameters also leads to the formation of classical 
bound electronic states in the MD simulations with no corresponding quantum counterpart. 

For high ion charges Z > 5, like N^+, Al^'^^ and Ca^°+ the PMFEX model is substantially improved and the 
agreement of the PMFEX-I- model with the MD simulations is excellent, as shown in Figs. [51151 and [TUl For a TCP 
with light ions, e.g. the case of Li'^+ ions in Fig. [H the PMFEX and PMFEX+ differ, but both somewhat deviate 
from the MD simulations as in the case of TCP. In this regime with moderate ionic charge state Z and for large 
coupling the renormalized Baranger-Mozer series should probably be truncated at higher orders. We have also shown 
that with increasing coupling the PMFEX scheme becomes invalid while the PMFEX+ approach still works and yields 
good or even almost perfect agreement with the MD, see Figs. [9l and fTOl 

In summary we note that the PMFEX approximation incorporates some key features of the MFD, in particular for 
large microfields E (small k values). The large fields are predominantly due to configurations with a single-plasma 
particle near the impurity ion. In this case the perturber-perturber correlations are less important and the PMFEX 
as an independent-particle model is accurate for both neutral and charged radiators. The behavior of the MFD at 
large E is closely related to the behavior of the generating function C {k) at small k and how much accurate the 
exact second moment is involved in the model. At small k the generating function behaves as C{k) ~ (k^/6)(£'^) 
(see Eq. ((HJ), where the second moment (i?^) is exactly involved in the PMFEX model. Furthermore the microfields 
in the large field regime are asymptotically Gaussian distributed and characterized by the second moment (i?^). As 
discussed above this is, however, not so for a neutral radiator and assuming bare Coulomb interactions between 
plasma particles where C (k) behaves as ~ k'^/^ at small k and the MFDs is similar to the Holtsmark distribution at 
large E. Conversely, the small microfields (large k values) are due to the additive effects of many particles at large 
distances and the PMFEX approximation is again accurate except in the case of the neutral point where some 
deviations from MD simulations may occur, see Figs. [3]15| and Ref. jl7| for other examples. The main uncertainty in 
PMFEX, therefore, remains the domain of the intermediate configurations with E 1 and k ~ 1 (in units of Eh 
and 1/Eh, respectively). Apparently, this domain is reduced by means of the renormalized (screened) fields £a (r) 
and the PMFEX is much more adequate compared to the standard Baranger-Mozer (second order) treatment. In 
PMFEX-I- model, on the other hand, imposing the renormalization procedure to the standard Baranger-Mozer series 
substantially improves their convergency, since ?/;a"^('^) tends to zero for the relevant configurations more rapidly 
than xi"\t^)- Finally, the PMFEX-f model may also be useful for a charged impurity ion in the cases where the 
PMFEX approximation deviates from MD. For the TCP, an example is the case of strongly coupled hydrogen with 
small regularization parameter 6 considered in Ref. p^ . 
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APPENDIX A: EXACT SECOND MOMENT FOR NEUTRAL RADIATOR 



In this Appendix we derive the exact second moment for the MFD in the TCP at the neutral point. We use the 
exact relation for the second moment derived in Rcf. 16] 

Ue{r)geR{r)dr- {r) g,n {r) dr\ , (Al) 

where Ua {r) = — [r^u'^j^, {r)]' with a — e,i and the pair interaction potentials UaR (r) are given by Eq. (jj]) with the 
regularization parameters 6a ji. 

The second moment {E'^) is ill-defined for the neutral-point distribution since Zn —> and gaR (r) — > 1 and the 
expression in the brackets in Eq. (|Aip vanishes in this case (this is true for the Coulomb potential as well as for any 
potential regularized at the origin) . To obtain the correct limit of Eq. (jAl|) at Zji we recall the definition of the 
correlation functions gaR {r) [16|. At vanishing this functions read 

gMr) - I ^ l^qaUaRir) +Y,<}pn0 J driUf^R{ri)[gai3{\r^ri\)^l]^+0{Zl). (A2) 



Here gap (r) are the equilibrium correlation functions. Note that the second term in the right-hand side of Eq. (|A2p 
is the excess potential energy of the TCP. In particular, using Eq. (|A2|) it is straightforward to calculate the effective 
field £a (r) for the neutral point. Insertion of this relation into Eq. (^5]) in the limit Zu — > and for Coulomb 
electron-electron, ion-ion interactions yield Eq. ((?7| with gaR (r) — 1 (see also Eq. 

Now we substitute Eq. (|A2p into Eq. (|A1|1 . In the limit of vanishing charge Zji this yields 

(A3) 



(E^) = 47rVn„ / El{r)r^dr 
Jo 



[gaa {r) - 1] (r) rdr -2 / [g^i (r) - 1] (r) rdr 



where Ea (r) = ~1aGFu'aR (^) is the single-particle electric field with = —1 and qi — Z , Up (r) = [rupn (r)]', and 

*a;3 {r)^l f [Up (\p-r\)-Upip + r)] UaR (p) pdp. (A4) 



The first term in Eq. (jA3[) is the averaged density of the electric microfield energy (self-energy), the second one is 
the density of the excess electric energy which appears due to the correlations between plasma particles. Assuming 
regularized interaction potential Eq. (jl]) we obtain Up (r) ~ e^'^^^'^'^ /SpR, and 

*e.(r) = l--e-^/^"^~- ^— , (A5) 

^ ' 2 2(1- A) 2(1 + A) ' ^ ' 



Here A = Siu/^eR- In this case Eq. (|A3|) yields 

— + — 



(E^) ^ 2nn^el ( ^ + ^ ) (A7) 

^ei {r) rdr 



{poo POO 
J2 J [5aa (r) - 1] (0 rdr ~ 2 j [ge^ (r) - 1] 



where 'i/aa (r) and 'i'^i (r) are now given by Eqs. (|A5p and (|A6p . The first term in Eq. (|A7p is precisely the second 
moment of the Hohsmark distribution obtained in Ref. [iB^ for the regularized interactions which is independent on 
Zfl. The second term arises due to correlations between particles. For the Coulomb interaction between plasma 
particles and the radiator with SeB. — or SiR ~ the first term in Eq. (jA7p diverges and the second moment of the 
MFD does not exist in this case. 
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